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The measurements of WMAP on the rehc density of dark matter strongly constrain 
supersymmetric models. In mSUGRA where the neutralino LSP is mostly a bino only 
rather fine-tuned models survive. On the other hand the relic density upper limit can be 
, easily satisfied in models with a Higgsino or wino LSP. 
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1 Introduction 



O 

I ' Despite the lack of direct experimental evidence for supersymmetry, the mini- 

Oh! mal supersymmetric standard model (MSSM) remains one of the most attractive 

extensions of the standard model. One of the nice features of the MSSM with R- 
parity conservation is that it provides a natural cold dark matter candidate, the 
lightest supersymmetric particle (LSP). Recently, the WMAP satellite has mea- 
sured precisely the relic density of cold dark matter, Qh^ = .1126^'Qi[g]^ (at 2a) 
$H ' This measurement severely constrains the parameter space of supersymmetric 

models. This is particularly true in mSUGRA models where most of the scenarios 
compatible with WMAP require a careful adjustment of parameters Indeed 
over most of the mSUGRA parameter space, the LSP is almost purely bino. As 
such, the LSP annihilates into fermions through the t-channel exchange of a right- 
handed sfermion. This process is not efficient enough to satisfy the tight WMAP 
upper limit on the relic density of dark matter. Then one must appeal to specific 
mechanisms such as rapid annihilation of neutralinos via s-channel Higgs exchange 
or coannihilation of neutralinos with other sfermions 0] to bring down the relic 
density in the desired range. Another possibility is that the LSP has a significant 
Higgsino component, then the couplings of neutralinos to the Higgs or Z is enhanced 
and annihilation into gauge bosons or fermion pairs becomes very efficient. One of 
these conditions is satisfied only for narrow strips in the parameter space of the 
mSUGRA model One basically finds three allowed region: the coannihilation 
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region which occurs usuaUy when m^^ « m^o ; the Higgs funnel region where anni- 
hilation via a heavy pseudoscalar Higgs into fermion pairs dominate at large tan /3; 
the focus point region where at large values of Mq one finds a neutralino with some 
Higgsino component. Note however that the location of both the focus point region 
and the Higgs funnel depend strongly on the standard model input parameters, in 
particular the top quark mass 5 , 6^ . 

The careful tuning of parameters necessary to comply with the upper limit 
on the relic density from WMAP within mSUGRA is however not generic of all 
MSSM models. Models with either a Higgsino or wino LSP generally satisfy easily 
the WMAP upper limit because of the large cross section for the annihilation of 
neutralinos into gauge bosons or fermion pairs T . Such scenarios can be found for 
example in SUGRA models with non- universal gaugino masses jS| El ■ Furthermore, 
in non- universal models, the relation between the heavy Higgs mass and the neu- 
tralino mass deviates from the mSUGRA prediction and one finds that annihilation 
through a heavy Higgs exchange can take place even at low tan/3. 

In order to assess the potential of future colliders to discover supersymmetry 
it is therefore important to consider the different categories of models and exam- 
ine the impact of the relic density constraint. Here we review the constraints on 
mSUGRA models showing also the impact of the top quark mass before discussing 
non-universal SUGRA models. Other constraints including direct searches as well 
as 6 ^ S7 are also briefly mentionned. 

2 Relic density of dark matter 

The computation of the relic density of dark matter amounts to solving the 
evolution equation for the density of the LSP which depends on the effective anni- 
hilation cross section [H| . Coannihilation processes involving SUSY particles slightly 
heavier than the LSP can contribute significantly to the effective cross-section. How- 
ever, these contributions are suppressed by a Boltzmann factor oc exp~^^/'^f where 
AM is the mass difference between the NLSP and the LSP and Tf the decoupling 
temperature. Here we will present results computed with the code micrOMEGAsp. 
This code includes all annihilation and coannihilation channels and takes into ac- 
count loop corrections in masses as well as in vertices. Loop corrections for the 
Higgs widths (including the Amf, correction) are also included. 

Within the context of a model defined at the GUT scale such as mSUGRA 
or non-universal SUGRA, one crucial step is the evaluation of the physical spec- 
trum corresponding to a given model. This is based on the renormalization group 
equations that specify the evolution of parameters from the high scale to the low 
scale as well as on the inclusion of radiative corrections to the physical masses of 
sparticles. The theoretical uncertainties in the prediction of the spectrum can se- 
riously affect the prediction of observables at low energies and in particular of the 
relic density Here we will mainly present results obtained with the spectrum 
calculator SOFTSUSYl.8.6 [Hlj interfaced with micrOMEGAs 1.3 through the SUSY 
Les Houches Accord [TT] . 
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3 Constraints on mSUGRA models 



We first remark that within mSUGRA the so-called bulk region (at low M1/2 
and Mo) that used to be one of the favoured regions has considerably shrunk after 
the much tighter upper limit on the relic density from WMAP. This is because in 
the bulk region the LSP is almost a pure bino. The annihilation of binos in fermion 
pairs is not efficient enough to satisfy the upper limit on the relic density of dark 
matter. One needs some additional contribution from the coannihilation or Higgs 
exchange channels. The important contribution of the latter at large tan/3 explains 
why the region allowed at low Mq — M1/2 is more important, see Fig.Q] Note that 
the low Mq — Ml/2 region is also constrained from the Higgs mass limit at least for 
intermediate tan/3 and by the measurement of 6 — > 57 when tan/3 > 35. We are 
then left mainly with three regions: coannihilation, focus point and Higgs funnel. 
Results for the case ~ 175 GeV are presented in each of these regions. 
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Fig. 1. Allowed regions in the Mo — M1/2 plane for a) tan/3 = 10, ^0 = 0, ^ > 0, 
mt = 175 GeV. Contours for for Qh'^ = .129(.094) (black/pmk dots), Mh = 111, 114 GeV 
(full green), Mh = 111 GeV,mt = 179 GeV (dash-dot), b) Allowed regions /or tan/3 = 50, 
Ao = Q, ^ > Q, mt = 175 GeV. The contour for b ^ s'y — 2.25 x 10"* is also shown 
(dash-dot). The LHC reach for L = 100/6~^ is reproduced from ll^Hl 

In mSUGRA, the coannihilation region designates the region at low AIq, where 
the lightest fi is the NLSP. In addition to neutralino annihilation into fermion pairs, 
the main coannihilation channels are Xi"^! ~* ^iid fifi — > rf. The relic density 
is very sensitive to the NLSP-LSP mass difference. The proper mass degeneracy for 
a relic density consistent with WMAP corresponds to a very narrow strip in the 
Mq — Ml/2 plane. In Fig.^ the allowed region correspond to the area between the 
Qh"^ = .094 and nh^ = .129 contours. Typically a mass difference AM^^o « 10 GeV 
for m^o « 
neutralinos, 



100 GeV is required while near degeneracy is necessary for heavier 
— 400 GeV. The mass difference between the LSP and the sleptons 



of the first two generations can be larger. Eventually, as the LSP mass increases the 
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coannihilation cross sections also become too small and one obtains an upper limit 
on the LSP mass. We find for tan/3 = 10 that m^o,mf^ < 448 GeV,mq < 2 TeV 
and rrig < 3 TeV. All this region should be within reach of the LHC . Within 
mSUGRA, one can also find models where coannihilation with the i-quark takes 
placets]. This occurs only for a large value of ^o- Although squarks are rather 
light, the main decay channel is into t —> cxi making it difficult for searches at 
LHC. 

The Higgs annihilation region includes both a light Higgs annihilation at low 
Ml/2 and a heavy Higgs annihilation at large tan/3. The former occurs only for 
the LSP mass very close to since the light Higgs has a very narrow width. 

The latter appears at moderate Mq — Mi /2 values and benefits from the enhanced 
couplings of the heavy pseudoscalar Higgs to 6-quarks. The main annihilation 
channels are XiXi ~^ f^b, tt corresponding to the preferred decay channels of the 
heavy Higgs. Due to the heavy Higgs annihilation channel, neutralinos as heavy as 
650 GeV can give reasonable values for the relic density. The only condition is that 
Ma — 2to^o < 2rA. This type of model could be very difficult to hunt at the LHC as 
it features a rather heavy spectrum. Indeed at the tip of the Higgs funnel, one finds 
all squarks and the gluino in the 2.5 — 3.2 TeV range. Note that the position of the 
Higgs funnel is sensitive to the standard model input value used in the spectrum 
calculator codes, in particular mb{m\j) 14 as well as, to a lesser extent, the top 
quark mass In evaluating the relic density in the Higgs funnel, the crucial 

parameters are M^ and m^o . Then the LHC can help constrain the relic density 
in this scenario with an accurate measurement of Ma if the A — > IJ.'^ i-J.^ channel 
can be used [T7j . 

The last allowed region is the focus point region, this region is found at high 
values of Mq where the value of n drops rapidly. This occurs near the electroweak 
symmetry breaking border. When n « Mi,M2, the LSP has a significant Higgsino 
fraction. The main annihilation channels are XiXi ~* W~^W~ , ZZ,tt. Coannihila- 
tion channels with heavier neutralinos and/or charginos are also important. The 
position of the focus point region is very sensitive to the value of the top quark 
mass that enters the RGE and also differs for different RGE codes, El ^| see 
Fig. Eb- For the case of SoftSUSY used here, rather low values of the top quark 
mass are necessary to reach the region where n drops rapidly when tan/3 = 10. 
For mt = 175 GeV, one finds a cosmologically allowed region at large Mq only for 
large tan /3 as displayed in Fig. QJi. Values consistent with WMAP for flh^ require 
a large XiXi^ coupling, that is not so large a value for fi. This implies that the 
chargino and lightest neutralinos should be rather light. For example, the chargino 
is confined to < 375 GeV for tan/3 = 50. Although the chargino and the LSP 
can be nearly degenerate this occurs mainly when the relic density falls below the 
WMAP range, for flh"^ « .1 typical mass differences are rather around 50 GeV. 
Note that models in the focus point region feature a very heavy sfermion sector. 
The LHC will have little opportunity to discover the squarks, only the gluino could 
be accessible. At the same time the rather light gaugino sector could be probed at 
a linear collider, extending the reach of the LHC ,12, . 
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Fig. 2. a) Allowed regions in the Mo — M1/2 plane for tan/3 = 50, Ao — 0, fi > and 
mt = 179 GeV. Same labels as Fig.l. h) Q,h^ vs mt for tan f3 — 10, Aq — 0, h > and for 
SoftSUSYl.8.6, Suspect2.3, Spheno2.2.0 and Isajet7.69 [7^ . 



3.1 Theoretical uncertainties: standard model parameters 

The values of the standard model input parameters play an important role in 
determining the SUSY spectrum in mSUGRA j^. Therefore they have an impact 
on the relic density constraint. This is true especially in the focus point and the 
heavy Higgs funnel regions. For one the heavy Higgs masses, which often provide 
the only efficient annihilation mechanism also depend on the top Yukawa as well 
as on mb(mf,) and as [51ll4|. More importantly the solution of the renormalization 
group equation for the ^ parameter, which drives the Higgsino nature of the LSP, 
is extremely sensitive to the top Yukawa 6 . 

A heavier top quark, say mt = 179 GeV, means a shift of the focus point region 
towards higher values of Mq in fact beyond Mq = 5 TeV even for tan (3 = 50. 
For large values of Mq, it is difficult to obtain a converging solution, although the 
authors of Ref. ^Sl find the focus point region all the way up to Mq = 10 TeV. The 
Higgs funnel moves closer to the coannihilation region when mt = 179 GeV. Indeed 
the heavy Higgs masses are shifted upwards with the increase of the top Yukawa 
making it increasingly difficult to have m^o « A! a/2. Note that increasing the bot- 
tom quark mass has the opposite effect. Finally the relaxed bound from the Higgs 
mass means that one recovers part of the region where light neutralinos/charginos 
and sleptons can be found. 



4 Beyond mSUGRA 

The nature of the LSP strongly influences the value of the relic density of dark 
matter so constraints on SUSY models can differ markedly from the mSUGRA case 
in models with a Higgsino or wino LSP. This can be achieved by relaxing the gaugino 
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universality condition. For one, the wino content of the LSP can be significantly 
increased in models where M1/M2 > 1 at the GUT scale. Second the value of fi can 
be much smaller than predicted in mSUGRA thus increasing the Higgsino fraction 
of the LSP. This occurs for example for models with M3/M2 < 1 at the GUT scale. 
Below we give explicit examples of these two types of models. Basically one finds the 
same mechanisms as in mSUGRA for getting sufficient annihilation/coannihilation 
of neutralinos but the allowed region in the parameter space of the MSSM shifts 
significantly and no longer requires a fine tuning of parameters especially if one is 
willing to allow models which are below the WMAP range. The direct constraints 
are of course also modified. 

4.1 Ml > M2 

Imposing Mi > M2 at the GUT scale lead to Mi > M2 at the weak scale and 
a LSP with a higher wino content. Then annihilation of neutralinos into W pairs 
becomes dominant. Since this process is much more efficient than the annihilation 
into ferniion pairs the relic density often falls below the WMAP range jZj . Consider 
for example a model with Mi = I.8M2. This value corresponds to Mi « M2 
at the electroweak scale. Then the LSP is a mixed bino/wino state. The allowed 
region in the Mq — M1/2 plane differs significantly from the mSUGRA case 0. 
The annihilation into gauge bosons occurs and has a large enough rate even when 
the fermion annihilation channels is small because the sfermions are very heavy. 
For example in Fig. O^i, the whole region to the left of the Vlh^ = .094 contour 
has a relic density below the WMAP range. Agreement with WMAP then implies 
a rather heavy LSP, for example m^o ~ 600 GeV for tan/? = 10. The tail in 
Fig. corresponds to a region where coannihilation with the chargino becomes 
significant, this coannihilation proceeds near the heavy Higgs pole. The constraints 
on this type of models from the Higgs mass as well as from precision measurements 
are not significantly different than in the universal model. 

As concerns collider searches, the spectrum of coloured sparticles is not dramat- 
ically different then in mSUGRA. However since the WMAP allowed regions are 
completely different and in particular include a Higgs funnel at small tan /3, models 
with very heavy squarks are perfectly acceptable. For example, for Mq = M1/2 = 
2 TeV and tan/? = 10, all squarks are in the 3-4 TeV range. These are beyond the 
reach of LHC. Furthermore, the LSP tends to be much heavier than in mSUGRA, 
more importantly, the small mass difference between the chargino and the neu- 
tralino can effect searches. Had we increased further the ratio M1/M2 we would 
have found a LSP with an even larger wino component. This would increase the 
values for the LSP masses compatible with WMAP. For example if M1/M2 = 2.5 
the LSP as well as the whole supersymmetric spectrum lies above 1.6 TeV unless 
there is an additional dark matter component. 

4.2 Ms < M2, the case A/3 = M2/2 

Decreasing the value of M3 while keeping M2 = Mi at the GUT scale also means 
that the relic density constraint can be satisfied more easily. First, one finds smaller 



A6 



Czech. J. Phys. 54 (2004) 



Relic density of dark matter in mSUGRA 



> 

o 



' 




tan/?=10 


_ V 










\ & 












1 1 i'! 1 : 1 ! 







5000 












4500 












4000 












3500 






3000 






2500 






2000 






1500 






1000 






500 






tan/?=10 



^2^ 



09* 



200 400 600 800 1000 1200 1400 1600 1800 2000 500 1000 1500 2000 2500 3000 

M,/2(GeV) M,/2(GeV) 

Fig. 3. Allowed regions in the Mq — M1/2 plane for tan/3 — 10, Aq — 0, jj. > and 
mt = 175 GeV a) Mi = I.8A/2 b) M3 = M2/2, Contours for Mn = 111, 114 GeF (green) 
and for Q.h^ = .129(.094) (black/pmk dots). 



values for which means a LSP with a significant Higgsino component and more 
efScient annihilation channels without requiring high tan/3 Second, one finds 

significant areas of parameter space where neutralino annihilation proceeds through 
s-channel Higgs exchange even for moderate values of tan f3. For example in the re- 
gion of a f NLSP, both stau-coannihilation and annihilation into fermion pairs via 
Higgs exchange contribute to the effective annihilation cross section. As a result, 
when the f — xl are nearly degenerate the relic density is already much below 
the WMAP range unless the LSP is very heavy. In this model, we also find rapid 
annnihilation trough a Higgs resonance but this time for the coannihilation pro- 
cesses XiX2>XiX3 ~^ ff or XiXi" ~^ tb via charged Higgs exchange. The onset of 
the charged Higgs exchange contribution is easily visible as a kink in the contours of 
constant relic density in Fig. |2|3. For large tan/3 the relic density constraint is sat- 
isfied in an even larger portion of the parameter space because of the more efficient 
annihilation/coannihilation through heavy Higgses as well as a large Higgsino-LSP 
region. Note that in these models both the Higgs mass and the 6 ^ 57 constraints 
become more severe as the squark mass are lower than in the universal case for a 
fixed value of M1/2 = M2\gut- 

A significant fraction of the models with M3 < M2 predict squarks within 
the range accessible by the LHC. However the region compatible with WMAP 
is not completely accessible even for intermediate values of tan /3. Because of the 
annihilation/coannihilation through the heavy Higgses, one can find reasonable 
values for the relic density even for very heavy neutralinos, hence a very heavy 
supersymmetric spectrum. For example for Mq — 4 TeV, = 3 TeV, the 

spectrum consists of heavy squarks, ruq = 3 — 4 TeV, neutralinos and charginos, 
m^o « « 820 GeV and also heavy Higgses mn/A ~ 1-5 TeV, a difficult task 
for discovering supersymmetry at colliders. The only new particle that could be 
reached at a collider would be the light Higgs. 
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5 Conclusion 

The latest measurements of WMAP on the rehc density of dark matter put 
strong constraints on supersymmetric models. In particular in mSUGRA very spe- 
cific relation among sparticle masses and parameters should be satisfied: neutralino 
nearly degenerate with a sfermion, neutralino mass near half the mass of a Higgs or 
neutralino with a significant Higgsino component. In mSUGRA most of the param- 
eter space would be covered at LHC except for the tip of the Higgs funnel region and 
part of the Higgsino LSP region. In non-universal SUGRA models, compatibility 
with the WMAP measurements imposes similar conditions on the SUSY spectrum 
than in mSUGRA. The main new possibility is for a wino LSP. However the relaxed 
constraints mean that a significant part of the allowed parameter space is beyond 
the reach of the LHC. 
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